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Abstract: 

Results of an investigation of the 0(4) spin model at finite temperature using 
anisotropic lattices are presented. In both the large N approximation and 
numerical simulations using the Wolff cluster algorithm we find that the 
ratio of the symmetry restoration temperature T$r to the Higgs mass mn is 
independent of the anisotropy £. From the numerical simulations we obtain 
a lower bound of Tg R /m H — 0.58 ±0.02 at a value for the Higgs mass mna s ~ 
0.5, which is lowered further by about 10% at mnfl s ~ 1. Requiring certain 
timelike correlation functions to coincide with their spacelike counterparts, 
quantum and scaling corrections to the anisotropy are determined and are 
found to be small, i.e., the anisotropy is found to be close to the ratio of 
spacelike and timelike lattice spacings. 
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1 Introduction 



The fate of a spontaneously broken gauge theory at finite temperatures of 
the order of the symmetry breaking scale has attracted attention for a con- 
siderable period now. Such investigations are of importance to the physics 
of the very early universe. Two prime examples are the inflationary uni- 
verse and the generation of the baryon asymmetry. It has been argued [|TJ 
that all the baryon asymmetry generated at the GUT scale is washed out 
by non-perturbative effects near the electroweak phase transition. Whether 
any extra mechanism exists to create a fresh baryon asymmetry || near this 
phase transition remains unclear. Although symmetry restoring phase tran- 
sitions in spontaneously broken gauge theories are crucial for these areas, our 
knowledge about them comes chiefly from perturbation theory || which can 
be expected to be rather inadequate for dealing with the anticipated presence 
of certain intrinsic non-perturbative effects near such phase transitions [|J. 

Motivated by the desire to learn more about the non-perturbative aspects 
of the symmetry restoration transition, exploratory lattice investigations of 
SU(2) Higgs-gauge models and nonlinear O(N) models at finite temperature 
have been made |5], ||, [7], || [!|. These models are expected to be trivial, giving 
rise to an upper bound on the Higgs mass of about 650 GeV at values of the 
(lattice) Higgs mass close to the lattice cutoff m-Ra s ~ 1. Correspondingly, at 
finite temperature one expects a lower bound on the symmetry restoration 
temperature Tsr in units of the inverse Higgs mass r/% . For the standard 
model the gauge couplings at the weak symmetry breaking scale and the 
Yukawa couplings, with the possible exception of that of the top quark, are 
small. Neglecting them as a first approximation, one arrives at an 0(4) 
symmetric scalar model. As in the case of the bound on the Higgs mass, one 
can then hope that the lower bound on Tg R /m H can be obtained by studying 
the 0(4) model, rather than the more complicated SU(2) Fermion-Higgs 
model. 

In this note we study the 0(4) model at finite temperature using nu- 
merical simulation of the euclidian path integral on lattices with anisotropic 
spacings in time (temperature) and spatial directions. We also compare our 
results with analytic results obtained in the lowest order large N expansion. 
Anisotropic couplings allow, at least in principle, a continuous tuning of the 
temperature, while the Higgs mass in units of the lattice spacing can stay 
fixed at values mna s ~ 1. Consequently a study of temperature effects of the 
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theory at a correlation length of the Higgs particle of order unity becomes 
feasible without changing the lattice size or losing the resolution in the tem- 
perature direction and the relevant information about the lower bound on 
T SR / m H can be extracted. In addition, anisotropic lattices allow us to distin- 
guish the finite temperature effects, which in the euclidian formulation that 
we employ could be regarded as a special type of finite size effects, from other 
finite size effects since the finite temperature effects have to be independent 
of the anisotropy in the scaling region. 

The plan of this paper is as follows: In the next section we define the 
model and give details of our methods to study it in the large N limit and, 
for iV = 4, using numerical simulations. The procedure to obtain the ratio 
Tsr/ttih is described here. Section 3 is devoted to the discussion of our re- 
sults and conclusions are presented in the final section. Some of our results 
have already been presented in a preliminary form in fTO . 



2 The anisotropic O(N) model 

The anisotropic O(N) symmetric spin model on lattices with spatial exten- 
sion N s and temporal extension N t is defined by the action 

s = -iV/% £ s x ■ s x+6 + i £ S x ■ S x+] ) (1) 

x 7 x,j 

or alternatively 

S = -2k( 7 E Sx ■ S x+6 + i £ S x ■ S x+] ). (2) 
x 7 ay- 

Here the spins S x are unit vectors in O(N), 7 is the anisotropy coupling 
and (3 or k denote the hopping parameter. Isotropic lattices are defined by 
7 = 1. For the study of the large N limit we take the first form of the 
action, eq. ([!]), keeping f3 finite, while for the 0(4) model we use the more 
conventional second form, eq. (0). 

Denoting the lattice spacing in spatial directions a s and in the temporal 
direction a t , the anisotropy parameter £ is the ratio of spacelike to timelike 
lattice spacings 

£ = -• (3) 

at 
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In the naive continuum limit and for noninteracting theories £ = 7. How- 



ever, quantum and scaling corrections can modify this relation |l l|]. For a 
given anisotropy coupling 7, £ can be determined by the requirement that 
physics, e.g., the fall-off of correlation functions, is the same in the time and 
spatial directions. As the 0(4) model is weakly interacting we expect only 
a small renormalization of £ with respect to the bare coupling 7. We also 
expect at the critical point £ = 7, as the renormalized coupling of the 0(4) 
model vanishes there. The anisotropy £ is easily calculable in the large N 
limit for the symmetric phase of the model. There we found that the relevant 
contributions to £ are of order 0(a 2 ), i.e., a scaling violation effect. The same 
conclusion can be inferred for the broken phase. We also looked at contribu- 
tions to £ in renormalized perturbation theory of the A0 4 model. Up to two 
loops we again found 0(a 2 ) effects, in contrast to theories involving gauge 
fields where quantum corrections of order 0(g 2 ) occur |TT[. We conjecture 



that for our model the difference of £ from 7 is order 0(a ; ) to all orders in 
perturbation theory. 

On the anisotropic lattice the physical 3-dimensional spatial volume and 
the temperature are respectively given by V3 = N^a^ and T = 1/N t a t = 
£/N t a s . It is therefore possible to vary £ and N t simultaneously at fixed 
ratio £/Nt, without changing the temperature in units a' 1 , or changing the 
spatial volume. This amounts effectively to a change in resolution in the 
time direction: at is changed while N t at is kept fixed. At least in the scaling 
region physical results should then be independent of the anisotropy £. A 
verification of this property will provide a valuable consistency check to our 
analysis. 

Earlier studies || |7| of the symmetry restoration phase transition in the 
0(4) symmetric spin model on isotropic lattices revealed that it was only 
possible to determine the symmetry restoration temperature Tsr for values 
of the Higgs mass which barely exceeded a value of m#a s ~ 0.4 on reason- 
able lattice sizes. Furthermore, increasing the Higgs mass in units of the 
lattice spacing a s , one expects a logarithmically slow decrease of T S R/m H , 
driving numerical simulations on isotropic lattices to larger temperatures and 
smaller Nt values, therefore loosing resolution in the time direction. However 
choosing the anisotropy coupling 7 > 1 it is possible to explore the model 
at values m#a s ~ 1 and at larger values of the temperature 1/N t a t without 
giving up a reasonable discretization in the time direction, i.e., in our case 
it was possible to simulate the region mna s — 1 on an N t = 4 lattice. In 
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this way it will be possible for the first time to explore regions of the theory 
where the Higgs mass takes values of the order of the cutoff and a numerical 
determination of the lower bound on T$r becomes feasible. 

In both the large N calculation and the numerical simulations our pro- 
cedure to investigate finite temperature effects consists of two steps. First 
we determine, at given value of the anisotropy coupling 7, the critical cou- 
pling on x N t lattices. Studying the large N limit in leading order, (3 C is 
obtained by solving numerically the saddle point equation 

for N s — > 00, where D(p) is given by 

D{p) = 4 7 W(^ ) +4£ S in 2 (~ Pi ) , (5) 

with the momenta p^ given by p^ = 27m M /iV M , n M = 0, . . . , — 1, where 
Ao = N t and Nj = N s . The prime on the sum in eq. @ indicates that the 
zero mode p = is being left out. In Monte Carlo (MC) simulations the 
unique crossing point of the Binder cumulant gn = (M 4 )/(M 2 ) 2 for various 
volumes N% and at given values of the anisotropy coupling 7 yields k c (oo, N t ). 
Here M is the order parameter, defined by M = ((M a M a ) - 5 ), where M a is 
given by 

M a = —J— V S% (6) 

and a denotes the O(N) index. Alternatively, one may use the peak position 
of the susceptibility x = N t N^((M 2 ) - (M) 2 ), to define k c (N s , N t ). Using the 
critical exponents of the 0(4) model in three dimensions, k c (oo, N t ) can then 
be obtained using the finite size scaling theory. We employed both methods 
and checked that they yield consistent results. 

Secondly the Higgs mass and the renormalized field expectation value 
were then determined at zero temperature at the coupling k c (oo, A^) on x 
^Ng lattices. For the determination of the renormalized field expectation 
Vr in units of we proceed in case of our Monte Carlo simulation as 
follows: The dimensionless quantity v-&a t is given by an estimator for the 
field expectation value E, which is properly normalized by its corresponding 
wave function renormalization constant Z: v^a t = YijyfZ. Note here that 
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(4) 



neither quantity £ nor Z are fixed numbers in the theory. It is possible to 
redefine £ and Z by overall multiplicative factors, such that the physical 
quantity v^a t stays fixed. In our case we chose the expectation value of 
the mean field multiplied with a convenient factor (v / 2k/7) < M > as an 
estimator for the field expectation value E. The corresponding wave function 
renormalization constant can then be derived from the behavior of the 0(4) 
symmetric zero momentum correlation function 

G W = Jj^E < SS^ne, >, (7) 

which is defined in the temporal direction of the lattice. Using chiral per- 
turbation theory one finds for large values of n on a periodic symmetric box, 
that G(n) has the shape of a parabola. This is due to the presence of massless 
Goldstone bosons in the theory: 

AT 

G(n) = Z—(n + const. (8) 

Expressing the volume V in units of a t , V = £, 3 N 3 N t , the desired wave 
function renormalization constant Z can in principle be determined. In our 
actual data analysis we have also considered the contribution of the scalar 
particle to eq. (|7|); for a detailed description of the procedure see |T2| . 



For the determination of the Higgs mass we project the scalar fields 
individually in each configuration onto the direction of the mean field M a / 
| M | and we obtain a field operator which has a good overlap with the Higgs 
particle: 

qa yrot 

| M | 

The Higgs mass m^dt can then be extracted from the exponential decay of 
the zero (spatial) momentum correlation functions of the operator S a>x . 

Introducing 0(N) invariant correlation functions defined on the main axis 
of the lattice in time direction 

cM) = -^-Y.s a x s: +net (io) 

^ * s t X 

and in space direction 

Cs(n) = J ^ w Y^s:s: +nea (ii) 
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we demand invariance with respect to an interchange of the spatial and tem- 
poral directions. We match the correlation functions in temporal and spatial 
directions at equal distance n, by scaling the temporal direction by a factor 
£, which determines the anisotropy. As we shall see below, the difference of 
£ from 7 was found to be rather small, being of the order of at most 3% for 
all 7-values we studied. 

3 Results 

The numerical computations have been performed using the nonlocal Wolff 
cluster algorithm. The employed statistics were about 10 5 sweeps for each 
simulated lattice size and set of couplings. At finite temperature we simulated 
N t and 7 values as given in Table 1. In each case we performed simulations 
with N s = 18 and N s = 24 at few values of the hopping parameter k. We 
employed the spectral density method in order to determine the maximum 
of the susceptibility and the crossing point of the Binder cumulant. At zero 
temperature, with k = n c (oo,N t ), we performed simulations on 18 3 x 7I8 
lattices with 7 equal to the cited values. 



N t 


7 


n c (oo,N t ) 


k c (oo, N t ) 


6 


1. 


0.3060(3) 


0.314594 


4 


1. 


0.3103(3) 


0.320871 


6 


1.5 


0.3645(3) 


0.377673 


8 


2. 


0.3912(3) 


0.408333 


3 


1.5 


0.3913(3) 


0.415998 


4 


2. 


0.4171(3) 


0.446373 



Table 1: Critical hopping parameters at given N t and 7 for the N s — > 00 
limit. The third row denotes our result from numerical simulations while the 
last row (k) denotes results from the large N expansion. 

Fig. 1 exhibits our results for both g R and x 011 18 3 x 6 and 24 3 x 6 
lattices for 7 = 1.5. We used the spectral density method to obtain the 
smooth curves shown from our data, shown by crosses. Similar results have 
also been obtained for all other values of 7 and N t . In each case we obtained 
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k c (oo, N t ) by using both the crossing point of qr and the finite size scaling 
of the peak position of the susceptibility. Both estimates were always found 
to be consistent, although we preferred to use the former for determining 
m H . Table 1 contains our results for function of N t and 7 from the 

numerical simulations, along with the corresponding results from the large 
N expansion. One finds a sizable but less than ~ 7% difference between the 
two estimates, which is of the same order as the discrepancy observed by 
comparing the zero temperature critical hopping parameter from the large 
N expansion with high precision numerical simulations. 

Fig. 2 compares the spacelike correlation function C s (n) on an 18 3 x 36 
lattice at (ft, 7) = (0.3912,2.0) with the corresponding timelike correlation 
function C t (n/£) at scaled distance n/£. One sees that the two are in nice 
agreement with each other. Table 2 contains, together with other quantities, 
the measured anisotropy £. The deviations of £ from 7 are small, on the few 
percent level, which is in accord with the expectations near a gaussian fixed 
point. 
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1.00(1) 


0.280(3) 
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.593(39) 
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.2082(2) 
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.1682(37) 


1.485(64) 
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.3645 


1.51(2) 


0.285(5) 


0.433(09) 


0.583(30) 










.1913(1) 


0.96(03) 


.1110(20) 


1.500(76) 
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.3912 


2.05(2) 


0.212(1) 


0.436(06) 
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.1879(1) 
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.0832(15) 


1.501(91) 


18 


27 


1.5 


.3913 


1.51(3) 


0.615(3) 


0.934(18) 


0.541(07) 










.3838(1) 


0.97(05) 


.2292(65) 


1.454(52) 


18 


36 


2.0 


.4171 


2.05(5) 


0.457(5) 


0.937(25) 


0.547(12) 










.3716(1) 


0.97(06) 


.1719(61) 


1.453(66) 



Table 2: Main results from the numerical simulation of the finite temperature 
0(4) model on anisotropic lattices. 
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Figure 1: Results for thermodynamic quantities at 7 = 1.5. 
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Figure 2: C s (n) and C t {n/£) at («,t)= (0.3912,2.0) on a 18 3 x 36 lattice. 
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The Higgs mass m^at was then obtained from an exponential fit to the 
connected zero momentum correlation functions of the operator eq. (9). 
These values of the Higgs mass are listed together with m H a<i in Table 2. 
Using these results, the ratio Tsn/mji = 1/N t mna t shown in the table was 
obtained for various 7 and N t . As expected, we observe the ^-independence 
of the ratio at fixed values of mna s , demonstrating the internal consistency 
of our finite temperature formulation of the theory on anisotropic lattices. 
Considering fluctuations around the saddle point in the large N limit, one 
can obtain the Higgs mass me at /3 c (oo, N t ) at given 7. Fig. 3 shows these 
large N results for Tsr/7Bh. They are also seen to be almost independent of 
the anisotropy £. 

We have also collected in Table 2 our results for the expectation value of 
the mean field < M >, the wave function renormalization constant Z, urg^ 
and finally the ratio Tsr/ur, which have been determined from the Monte 
Carlo data by the methods described above. Once again the £ independence 
of the ratio at fixed mna s is nicely born out. Large N results for the ratio 
Tsr/vr are also shown in Fig. 3. The renormalized vacuum expectation value 
of the field in the large N calculation is given by t> R = N{(3 c {^L t ) — (3 c (oo)) 
and we have set N = 4. 1 Again, the ratio Ts R /u R is almost independent of 
the anisotropy £. The large N results, shown in Fig. 3, agree quite well with 
the numerical results at N = 4 of Table 2. 

A remark concerning the error determination for the quantities as cited 
in Table 2 and a comment on further possible systematic errors may be ap- 
propriate here. As can be noted, the ratios T SR /m H and T SR /t> R exhibit 
sizable errors, as compared to the relatively small and purely statistical er- 
rors quoted for all the other quantities. These errors are mainly caused by 
the uncertainty of the finite temperature critical k values (Table 1), which in 
turn lead to relatively large errors for the zero temperature values of m-nat 
and v-RCit to be used in the ratios. Also we have to expect zero temperature 
finite volume corrections to the quantities m R at and v^a t used to construct 
the ratios as quoted in Table 2. As we anticipate the finite volume corrections 
to the quantities Tsn/mn and Tg R /t> R , quoted in Table 2, to be significantly 
smaller on our lattices than the errors induced by the uncertainty of the 

1 The normalization of vr in (To) differs by a factor (= 2 for N = 4) from the one 
used here. This causes a difference of a factor 2 in the scale of Fig. 2 there as compared 
to Fig. 3 in this paper. 
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Figure 3: Large N results for T SR /v R and T SR /mu as a function of the 
anisotropy parameter 7. 



11 



critical points, we refrained from a detailed zero temperature finite size scal- 
ing analysis. Future simulations yielding more precise k c values will have to 
incorporate them. 

Our data for Tg R /mH as depicted in Table 2 decrease, as expected, very 
slowly as the Higgs mass mHO s in units of a s is increased. Thus, depending 
on the choice of value of the correlation length up to which an effective theory 
can be defined, one obtains a lower bound on the ratio T SR /m H . Just as in 
the case of the determination of the upper bound to the Higgs mass, it is 
expected that this lower bound saturates for the theory under study, i.e., 
the 0(4) model at infinite bare quartic coupling. From Table 2 we estimate 
this bound to be 0.58 ± 0.02 for a correlation length of ~ 2, which further 
decreases by about 10% for a value of mHO s ~ 1. Our data for Tsr/v-r 
show an approximate constant behavior as mna s is varied. The actual value 
is within the errors consistent with the value ^/2, which is the prediction 
of one-loop renormalized perturbation theory, though the data show some 
tendency to lie slightly above the perturbative value. 

It is interesting to compare our results for Tsn/mn with the one- loop 
result as obtained in renormalized perturbation theory in the 0(4) model. 
To this order the symmetry restoration temperature is given by [0] 

^ = (-)*, (12) 
m H g R 

where gn is the renormalized quartic coupling of the model. Using previous 
high precision numerical determinations of #r |13| as an input we draw in 
Fig. 4 our numerical results for Tsn/m-^ (crosses) together with the one-loop 
prediction as indicated by the curve and by the triangles. Here we observe 
sizable deviations when mna s takes values ~ 1, indicating that higher order 
corrections are large at finite temperatures in a region of the model where 



the scalar correlation length is close to 1; see also [p~4| . Including also results 
from the large N expansion in Figure 4, one also notices sizable deviations 
of the large N results from our data, though the overall trend is reproduced. 

4 Conclusions 

Using anisotropic lattices we have studied the finite temperature behavior 
of the 0(4) theory in regions of the parameter space where the correlation 
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Figure 4: TsR/mn as a function of mna s . The crosses denote our numerical 
results, circles and boxes come from the large N expansion at various values 
of N t with anisotropy parameter 7 = 1 (circles) and 7 = 2 (boxes), while the 
curve and triangles correspond to 1-loop renormalized perturbation theory, 
eq. (12). 
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length of the scalar particle is as low as ~ 1. Depending on the maximal 
value of muds one is willing to admit for a sensible definition of the effective 
theory, a lower bound on T sr /toh is derived. E.g., for a heavy Higgs particle 
which at a value of the cutoff mno s ~ 0.5 has a mass close to its triviality 
bound of about 650 GeV, we find Tsr. = 370 GeV. This value is close 
to the value predicted by renormalized perturbation theory Tsr = v^^weak 
with f wca k — 250 GeV and consistent with our finding that the ratio T SR /t> R 
follows the perturbative answer in the whole considered correlation length 
interval. However, at correlation length 1 we start finding large deviations 
from one-loop perturbation theory for the quantity Tsr/ttih. Qualitatively, 
the lowest order large N expansion seems to reproduce all the features of the 
Monte Carlo (MC) data well. Even quantitatively the results are consistent 
with the naive expectation that they should be accurate to 0(1/ N). In the 
large N expansion we were able to explore ^-independence of Tsr/ttih and 
Tsr/^r over larger ranges of £ and for more values of N t . This supports our 
belief that the early scaling evidence in the MC data even for N t = 3 and 4 
lattices is no fluke. But it would be interesting to check this by simulating 
the theory at more £ and N t values. 
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